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We investigate dynamical chiral symmetry breaking in the Coulomb gauge Hamiltonian QCD. 
Within the framework of the linked cluster expansion we extend the BCS ansatz for the vacuum 
and include correlation beyond the quark-antiquark paring. In particular we study the effects of 
the three-body correlations involving quark-antiquark and transverse gluons. The high momentum 
behavior of the resulting gap equation is discussed and numerical computation of the chiral symmetry 
breaking is presented. 
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I. INTRODUCTION 



Chiral symmetry plays a major role in constraining 
the spectrum of low energy QCD. At zero density it 
is spontaneously broken and the associated Goldstone 
bosons dominate the low energy, soft hadronic interac- 
tions. The quark-giuon interactions which in vacuum 
break chiral symmetry may in dense matter, e.g. in 
the interior of neutron stars, lead to other, novel phases 
of the quark gluon plasma [Q. The chiral properties 
of the QCD vacuum at zero temperature and density 
have been extensively studied in various approaches to 
soft QCD g |, I I |. In principle one could in- 
vestigate it using lattice gauge methods. However, ex- 
trapolations of lattice simulations to small quark masses 
iTiu.d << 50—100 MeV (chiral extrapolation) still present 
a major challenge. In approaches bases on a Dyson- 
Schwinger formulation of QCD, dynamical chiral sym- 
metry breaking can be studied by analyzing the be- 
havior of the quark propagator. Even though no sys- 
tematic truncation scheme of the Dyson series in QCD 
exists, and in a majority of studies model interactions 
are introduced, the approach gives a good description 
of the low energy phenomenology. In particular it en- 
ables to correctly predict many of the static properties 
of the low lying mesons and baryons, i.e. masses and 
charge moments, and simultaneously account for the dy- 
namical chiral symmetry breaking as measured by the 
vacuum expectation value of the scalar quark density, 

{^PtP) (250 MeV)3 0. This value follows from PCAC, 

Goldstone theorem and the current algebra which result 
in the Gell-Mann-Oakes-Renner (or Thouless theorem) 
relation, —2mq{^p^p) = f^m^. Here, to, ^ 5 — 10 MeV 
is the current light quark mass, renormalized at the 
hadronic scale, /^r = 93 MeV is the pion decay constant 
and m^r is the pion mass. Without explicit chiral sym- 



metry breaking niq — 0, the above relation cannot be 
used to determine {^ptp)■ However, as rUg no phase 
transition to a chirally symmetric state is expected, and 
therefore the - (200 - 250 MeV)^ should still be a good 
estimate of the condensate in the chiral limit. 

Spontaneous chiral symmetry breaking enables to put 
the constituent quark representation of hadrons in a firm 



theoretical ground. The bare quark states defined with 
respect to the perturbative vacuum are replaced by quasi- 
particle excitations of the chirally noninvariant ground 
state. Residual interactions correlate the quasiparticles 
to form composite hadrons in which each valence quasi- 
particle contributes kinetic energy of the order of a few 
hundred MeV. This is analogous to the constituent quark 
model representation of hadrons and therefore it might 
be possible to further constraint the quark model phe- 
nomenology from a first principle, QCD based analysis 
of dynamical chiral symmetry breaking. Since the quark 
model picture calls for a Fock space representation it is 
most natural to consider a canonical, time-independent 
formulation of QCD. Coulomb gauge QCD offers such a 
framework [|[ ^ |l^. In the Coulomb gauge the single 
particle spectrum contains only physical degrees of free- 
dom, i.e. two transverse gluon polarizations. As long as 
the gauge fields are restricted to the fundamental mod- 
ular region, with no Gribov copies, the Hamiltonian is 
positively defined, it leads to a continuous time evolu- 
tion, and it is amenable to a variational treatment. Fi- 
nally the Coulomb gauge formulation leads to a natu- 
ral realization of confinement. This arises because elim- 
ination of the non-physical degrees of freedom through 
the gauge choice, V • A = results in an effective, long 
ranged instantaneous interaction between color charges. 
This interaction is the analog of the Coulomb potential in 
QED. In QCD however, the colored Coulomb gluons can 
couple to transverse gluons leading to a Coulomb kernel 
which also depends on the dynamical gluon degrees of 
freedom. As shown in Ref. |ll[| summation of the domi- 
nant IR contributions to the vacuum expectation value of 
the Coulomb operator results in a potential between color 
charges which grows linearly at large distances in agree- 
ment with lattice calculations ||l^. In a self-consistent 
treatment the same potential modifies the single gluon 
spectral properties and leads to an effective mass for 
quasi-gluon excitations 0(500 — 800 GeV), which is also 
in agreement with recent lattice calculations. The ap- 
pearance of the gluon mass gap can be used to justify 
the implicit assumption of the quark model that mixing 
between valence quarks and Fock space sectors with ex- 
plicit gluonic excitations is small. We will return to this 
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point in Section III. 

The Coulomb gauge formulation provides a very nat- 
ural starting point for building the constituent represen- 
tation in accord with confinement and dynamical chiral 
symmetry breaking. However, as it was noticed some 
time ago in the Coulomb gauge the simple BCS treat- 
ment of the vacuum is not sufficient to generate the 
right amount of chiral symmetry breaking. In particu- 
lar if a pure linear potential is used, V{r) = br with 
b ^ 0.2 — 0.25 GeV^ as determined by lattice calcula- 
tions one typically obtains U'M)]^^'^ ~ 100 MeV i.e. too 
small by a factor of two ^ . The short range part of 
the Coulomb potential requires proper handling of UV di- 
vergences and renormalization and in most recent studies 
has been ignored. As will be shown later, it does signif- 
icantly enhance the condensate and we will argue that 
the missing contribution can be accounted for by three- 
particle correlations on top of the BCS-like, particle-hole 
vacuum. 

The paper is organized as follows. In Section II we 
briefly discuss the canonical Coulomb gauge formalism 
and the linked cluster expansion which enables to include 
multi-particle correlations into the many-body ground 
state. We will derive the resulting contributions to the 
mass gap including up to three-body correlations. The 
formalism is suitable for handling both zero and finite 
density system and in this paper we will focus on the 
former. In Section III we discuss the approximations, 
numerical results and possible sources of UV divergence 
and their renormalization. Our conclusions and outlook 
are given in Section IV. 



II. COULOMB GAUGE HAMILTONIAN AND 
THE LINKED CLUSTER EXPANSION 



QCD canonically quantized in a physical gauge, e.g. 
Coulomb gauge, results in a Hamiltonian that can be 
represented in a complete Fock space defined by a set of 
single particle orbitals. One possibility is to choose the 
single particle basis as eigenstates of the kinetic (nonin- 
teracting) part of the full Hamiltonian, 

Ho = H{g = 0)^ I dxV/(x)[-i- V + /3m]^^(x) 



dx TrH(x)2+Tr (V x A(x))-^ . (1) 



The vacuum, |0), of Hq is shown schematically in Fig. la. 
The singe-particle excitations at zero density correspond 
to adding gluons to the positive energy, parton-like lev- 
els and quark-antiquark paris by creating a particle-hole 
excitation around the zero-energy Fermi surface. These 
excitations have energies given by, eg(k) — e?(k) = 



canonical commutation relations for transverse fields, i.e. 

[H'^(x), A\y)\ = ~i5'^'5t{V)5^{^ ~ y), (2) 

where (5t(V) = / — V® V/V^. In terms of the single par- 
ticle creation and annihilation operators, the color triplet 
of quark fields (i = 1, 2, 3) is given by. 



A=±l/2 



dk 

(2^ 



[u(k,A)6(k,A,i) 



-t;(-k,A)d^(-k,A,i)] e**'-". 



(3) 

where u and v are solution of the free Dirac equation for 
a fermion with mass m. In the following we will restrict 
our discussion to chirally symmetric case i.e. from now 
on we will set m ~ Q. The giuon field is given by, 



1 



[a(k, A, a)e(k, A) 



V»(k) 
+at(-k,A,a)e*(-k, A)] e^'^'", 

(4) 

with a;'^(k) — eg(|k|). The unrenormalized Coulomb op- 
erator is given by, 

^2 



c 



J dxdyp° (x)i^,fc(x,y,A)p''(y), 



(5) 



where p°(x) = ■(/'^(x)T°i/'(x) + f'^^W'ix) ■ A'=(x) is the 
color charge density and the kernel K is given by. 



Kab(:x.,y,A) = (x, a| 



1 



V D 



1 



VD 



|y,^'>, (6) 



where D is the covariant derivative in the adjoint repre- 
sentation, and the (x, a\ ■ ■ ■ |y, b) matrix element is given 
by (x,a|D|y,6) = [S^'V^ + gr^''A^{x)] 5^jc-y), and 
(x, a|l/V^|y, 6) = — l/47r|x — y|. When He is normal 
ordered with respect to the perturbative vacuum, |0) one 
might expect that the mean field, Hartee-Fock correc- 
tions to the single particle energies could already gener- 
ate an effective mass. This is not the case. The vacuum 
is a color singlet and thus the direct contribution from 
He to a single fermion energy vanishes. Furthermore, 
chiral symmetry of the Hamiltonian and of the perturba- 
tive, |0) vacuum protects the exchange term from mass 
generation. The effective mass can only be obtained 
if quark-antiquark correlations are introduced into the 
ground state as shown schematically in Fig. lb. 

The full, unrenormalized Coulomb gauge Hamiltonian 
has the following structure jl0[ ^ ^ 



H = Ho + Hc + Vgg + V3g + Vi^ 



- |k| for quarks, antiquarks and glu- YLere ¥„„ is the quark-transverse gluon interaction, 



ons, respectively. The quark fields in Eq. satisfy the 
canonical anticommutation relations and the gluons fields 
are given by H = HT" and A = AT" and satisfy the 



^99 = 5 j dxip'^{x)a ■ AV'(x), 



(7) 



(8) 
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and V^g and V43 represent 3- and 4- transverse gluon 
couplings arising from the nonabelian part of the mag- 
netic field, B'^ = V X A*^ + gP^''A^{yi) x A'=(x). Fi- 
nally Hcorr contains terms which come from a commu- 
tator of the determinant of the Faddeev-Popov operator, 
= Det{\J ■ D) and the gluon canonical momentum H. 
The detailed analysis of this Hamiltonian, emergence of 
confinement and issues related to renormalization in the 
gluon sector were discussed in Ref . . 



FIG. 1: Schematic representation of the particle, n"*" = (6^6) 
(solid) and hole n~ = 1 — {d^d) (dashed) occupations in par- 
tonic (left) and BCS (right) ground state 



A. Linked cluster expansion 

Since the Fock space basis generated by the set of single 
particle creation operators, 5^, d\ is complete, the true 
ground state, of H can be written as, 



\n) = 



(99; 
12 



12 



123 



1234 



-f^l234 Ol"203"4 



|0). 



(9) 



Here Ji2...„ represent wave functions of n-body clus- 
ters in the vacuum, and 1, 2 • • • collectively denote quan- 
tum numbers of single particle orbitals. This expan- 
sion is however impractical since it does not differenti- 
ate between connected (linked) and disconnected con- 
tributions. For example, at the 2-quark-2-antiquark 
level there are disconnected contributions of the type, 
F[llf^ = F[f''F^^f\ I.e. part of the n-particle cluster 
contribution originates from products of smaller, m < n, 
TO-particle clusters. 



The essence of the linked cluster expansion is based 
on the observation that all multi-particle correlation in 
the ground state, including the disconnected ones can 
be accounted for by proper resummation of the linked 
clusters only. This is achieved by writing the full ground 
state as jl^ 

\n)=e-'\0), (10) 
with S having the expansion 



12 



0(9 
012 



'a\a2 + 



^llT^^ldWa + E '^123r^^l'^2^3'^4 



123 



1234 



(11) 



with the operators S including connected pieces only. 
Comparing Eq. (||) and Eq. (O) we find for example that. 



ones. The expansion coefficients, s[2\.„ can be deter- 
mined from the eigenvalue equation for \ 



F, 



im) 
12 



S- 



iqq) 

11 ! 



123 



q(qqa) 

"^123 ' 



p(9999) _ 0(9999) I }_ 0(99) 0(99) 
-'^^1234 —'-'1234 "T 2 12 34 ' 



(12) 



i.e. the general expansion of Eq. (||) is obtained with all 
disconnected contributions constrained by the connected 



e'^He-^lO) = En\0). (13) 
This equation projected onto the partonic Fock space ba- 

(n) 

sis leads to a set of equations for the amplitudes Si2...n 
and the ground state energy, Eq, 



J 



{0\e'^He ^\0)=En, {qi,q2, ■ ■ ■ qn,;qi' , , ■ ■ ■ Qn,; 9i" , 92" , ■ ■ ■ 9n,\e^ He '^{O) = 0, Uq, rig, Ug = 1,2, ■■■ . (14) 
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In a nonrelativistic many-body system the Hamiltonian 
is typically a polynomial in the field operators. Since 
each S*'"' contains only particle creation operators, the 
matrix elements of He~^ between an n-particle state 
and the free vacuum will involve only a finite number 
of terms arising from the expansion of the exponentials. 
For example in a typical case when H — Hq + V with 
Hq being a one body {e.g. kinetic) operator and V a 
two-body potential one has, 



H + [S, H] 



[S,[SAS,[S.H]]]]. (15) 



In this case an approximation to Eq. (p^, is fully speci- 
fied by a number of clusters retained in S. This is, how- 
ever, not the case for the relativistic system discussed 
here. The expansion of the Coulomb kernel leads to an 
infinite series of operators to all orders in the transverse 
gluon field. Thus an approximation to Eq. (|lj) consists 
of specifying which clusters are kept in the definition of 
S and of a truncation scheme in evaluation of matrix 
elements of e^He~^ . 

The truncation of S limits the number of quark- 
antiquark-gluon correlations build into the ansatz for the 
ground state. At first one might think that such a trunca- 
tion would be hard to justify since any hadronic state, in- 
cluding the vacuum should have a large (infinite) number 
of partons. However, the first two terms in S, S'^'^^) and 
S'(»s) change the single particle excitation spectrum and 
effectively replace the partonic basis by that of massive 
quasiparticles. This is known as the Thouless reparam- 
eterization and is equivalent to the BCS ansatz for 
the vacuum which contains two-body, quark-antiquark 
and gluon-gluon correlations. The BCS ansatz leads to 
the chiral gap, constituent mass for the quarks as well as 
effective mass for the transverse gluons. Iterative contri- 
butions of multiparticle states which determine the wave 
functions of larger clusters, 5"*^"^ n > 2 are therefore 
suppressed by the quasiparticle energy gap. This gap is 
0(400 — 600 MeV) for quark-antiquark excitations and 
0(500 MeV - 800 GeV) for a gluonic excitation. The 
former follows from the typical constituent quark mass 
and the later from the gluon spectrum in a presence of 
static color sources as calculated on the lattice [n2| and 



are consistent with explicit calculation using the BCS 
gluonic ansatz for the Hamiltonian The transfor- 

mation form the partonic to the quasiparticle basis, gen- 
erated by 5"^^^ proceeds as follows. The (unnormalized) 
quasiparticle, BCS vacuum \ flBCs) is defined as, 



irises) ^e-^'''|0). 



(16) 



ith, 



^(2) ^ ^(99) + ^(99) = J2 S[fb{4 



12 



(99) J J 



SO that 



A canonical transformation which maps the set of free 
particle operators b, 6^, d, S , a, onto a set of quasipar- 
ticle operators B, B\ D, D\a, is defined by 



(17) 



(18) 



Bi = 
Di = 
ai 



1 



S 



1 

^l+\Si99)\^ 
1 

y/l-\S(99)\' 



2 



(gg) 

12 



, ^l + \Si9g}\2 

c(99) 
'^21 



-Ml 



E 



s 



iaa) 
12 



yJl-\S(99)\ 



r4, (19) 



where \S(9g)\'^5i2 



S(gg)S(9g)^ 



and similarly for 



|5'(9s)| These quasiparticle operators satisfy the canoni- 
cal (anti)commutation relations, they annihilate the BCS 
ground state. 



Bi\ilBcs) = Di\nBcs) = cti\nBcs) 



0, 



(20) 



and generate a complete Fock space. The eigenvalue con- 
ditions for the vacuum, Eq (Q) can therefore be rewritten 
in the quasiparticle basis. 



J 



{nBcs\e^He-^\nBcs) 



BCsl^BCS 



= Er, 



{Qi 



}nQ',Ql',Q2':---QnQ',Gl",G2", ■ ■ ■ Gnol^ He \^BCs)—0 

nQ,nQ,nG = 1,2,---. 



(21) 



Here the operator S contains contributions from 3- 
quasiparticle cluster and higher. 



j:s[g^-^BlDlal 

123 



(22) 



The matrix elements «Sj2''-.n can be related to S'j2.''..n by 
replacing the free particle operators by the quasiparticle 
operators. From the structure of Eq. ( [l9| ) it follows that 
for given n the operators 5'"' are a linear combination 
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of 5'^*^ including i < n. Since Eq. (|l9|) defines a canon- 
ical transformation the two sets of equations, Eq 
and Eq. are equivalent and one can simply use the 
later i.e. work directly in the quasiparticle basis with- 
out referring to the partonic basis. As suggested by the 
quark model it is preferred to represent low energy QCD 
eigenstates in terms of quasiparticle, quark and gluon ex- 
citations. From now on we will work the matrix elements 
of S in the quasiparticle basis and for simplicity rename 
them as ^ 

As mentioned earlier, in QCD, with S = 5^"-' trun- 
cated at some maximal n, Eq. ( ^l|) still contain an infinite 
number of terms arising from the expansion of He^^ . 
Since this (infinite) series is related to the multi-gluon 
structure of the Coulomb operator, _ftr(x, y,A), it can 
be organized according to how each of the terms renor- 
malizes the 0-th order Coulomb potential, Kab{x, y, 0) = 
^a&/47r|x — y|. To illustrate this consider truncating S at 
n = 2. The Ihs. of the first equation in ( pl| ) reduces to 
the expectation value of H in the BCS vacuum, 

{nBcs\e^He-^\nBcs) = {ilBCs\H\nBCs) ■ (23) 

The lowest order (in the loop expansion) diagrams are 
shown on the left side of Fig. 2. The matrix element 
{^BCs\H\flBCs) defines an effective potential ^//(x — 
y), by 

Veff (x - y) s {ilBcs I K^ti^, y , A) | ^bcs) • 

(24) 

It is straightforward to identify diagrams which give the 
dominant contribution to Ve// in both, the IR (|x— y| >> 
1/Aqcd) and the UV (|x - y| « 1/Aqcd)- In the 
IR region these are given by diagrams which, at a given 
loop order contain the maximum number of soft poten- 
tial, if(x, y,0) fines; the UV region is dominated by 
loops with the smallest number of vertices. The series 
of ring and rainbow diagrams, shown in Fig. 2, accounts 
for the leading IR and UV contributions to K//, re- 
spectively. The approximation can be systematically im- 
proved by taking into account the subleading contribu- 
tions e.g. vertex renormalization jlTl H If larger clusters 




FIG. 2: A few lowest order contributions to Ve//. Tfie 
two right most diagrams represent diagrams in the ring and 
rainbow series. The dashed line stands for the 0-th order 
Coulomb potential g^Ko = (?^/47r|x — y|. 

in S are retained, the expansion of e^ He~^ generates op- 
erators that have nonvanishing matrix elements between 
the vacuum and states with an arbitrary large number 



of particles. This is because as long as S'^"-' contains 
a gluon operator an infinite number of commutators, 
[S'("),[S'("),[---[S'("\i/] •••]]] arc nonvanishing. Their 
contribution arise from contracting gluons from each S**^") 
with gluons from the Coulomb operator. For example a 
term in S which contains pure glue operators (no quark 
or antiquark) will contribute to any matrix element in 
Eq. ( |2l| ) with any number of particles (gluons). This 
is illustrated in Fig. 3 for S^^^\ It is clear, however, 
that this type of corrections have the effect of simply 
renormalizing Ve//, beyond the BCS-like contributions 
shown in Fig. 2. Since the operators S'*^"^ commute with 
each other one possibility is to consider the effects of 
the pure gluon operators first, generate the new effec- 
tive interaction and then introduce clusters which con- 
tain quark and antiquark operators. Since H is a finite 
order polynomial in the quark operators, each term in S 
containing only quark and antiquark operators will lead 
to a finite number of terms in a matrix element between 
the BCS vacuum and a multiparticle state with a fixed, 
n = ng + ng + no- 




FIG. 3: Example of contributions to 

[5'("',[5'<"),[---[S("', //]•■■]]], for S = 5''""' and H = He- 
The two diagrams contribute to a matrix element with 
no = 2 and no ~ 4 respectively. Here the dashed line repre- 
sents the Coulomb potential dressed by the BCS corrections, 
Eq. (^^, e.g. through resummation of the ring-rainbow 
series shown in Fig. 2. 

To summarize, the linked cluster expansion of the QCD 
ground state is much more complicated that in a typical 
nonrelativistic many-body problem. Nevertheless it can 
be used to systematically improve the BCS approxima- 
tion. It is important to notice, however, that even the 
BCS ground state already probes the nonabelian multi- 
gluon dynamics via {QBCs\K{x,y, A)\nBCs) ■ In BCS 
this leads to and effective interaction which is very close 
to the potential between color sources and when treated 
selfconsistcntly leads to a quasiparticle (constituent) rep- 
resentations. 

B. QQG contribution to the quark mass gap 

In the following we will concentrate on the dynami- 
cal chiral symmetry breaking and therefore consider vac- 
uum properties in quark sector. As mentioned earlier the 
BCS mechanism of quark-antiquark pairing seems to be 
insufficient to account for the full dynamical symmetry 
breaking. We will discuss this point quantitatively in the 
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FIG. 4: Operators from Eq. |27| which contribute to matrix 
elements (nq, ng, ng | ■ ■ ■ |ri_BCs) for tlq < 1, Uq < 1 and 
no < 2 As in Fig. 3, the potential (dashed) line is the Veff 
of Eq. (|2^). The matrix element, S^^QG) corresponds to the 
rightmost vertices. 



following section. Our interest here is in extending the 
BCS approximation by including the effects of the next 
to leading (beyond BCS) order in the cluster expansion 
i.e. the 3-particle cluster contribution to the vacuum. 
We will therefore study, 

^_,^(QQG)^^5(QQG)^t^t4. (25) 

123 

The quark gap equation follows from, 

(Q,Q,|e^"^^^'i/e-^''"'^'|17,c5)=0. (26) 

This equation determines single particle orbitals and 
therefore it also gives the quasiparticle spectrum via, 
£1^12 = (Qi|^^|Q2)- There is a finite number of terms 
contributing to, Eq. ( p6[ ) 



= {QiQ2\H + [^(Q^G)^ ^] + ^[^(WG)^ [S^QQ'^Kh]] + |[5(WG)^ [S^QQG)^ [S^^'^'^K H]]\nBcs). (27) 



r 



The series is finite because starting at 0{ 



SiQQG) 



) commutators will produce operators which 
have at least 2-quark and 2-antiquark creation operators 
and these vanish between {QQ\ and | rises)- Some of the 
contributions to Eqs. (|2^) and (^) are shown in Fig. 4. 



In order to solve Eq. ( |27| ) and determine the single par- 
ticle basis, it is necessary to first solve for the amplitude 
5(Q0g) xhis amplitude can be obtained by projecting 

q(QQG) o(QQG) 



onto the three particle cluster. 



3!' 



BCS) 



0, 



(28) 



which also contains a finite number of terms. The two 



equations Eq. (27) and Eq. ( pSD form a set of coupled non- 
linear, integral equations for the amplitude S'('3'3G) g^j^j^ 
the single particle orbitals (or the BCS angle, Eq. (p^). 
In this paper we will simplify t^hese equations by lineariz- 
ing them with respect to 5''^'3QG)^ Eq.(|2|) then yields. 



I 



En — 



{n\H\nBcs) (29) 



Here \n) is the set of cigenstates of H in the three particle 
QQG subspace. 



H\n) = {E^^EnscsM 



(30) 



The contribution from S'^'^QG) ^j^g quark gap in 
Eq. ( p7| ) is then given by, 

5i2,n5mg = {QiQ2\[S^'^'^''\H]\nBcs) 



Y.{QiQ2\H\ 



1 



En E[2gf-,g 



{n\H\nBcs)- 
(31) 



Here, (5i2,n symbolizes the product of all (5-functions 
which restrict the quantum numbers of \Q1Q2) to be 
same as of the vacuum. With inclusion of Snig the gap 
equation can be written as. 



= SniQ + 5mc + Srug, 
where the BCS part given by 



(32) 



h2,ii[5mQ 6mc] = {QiQ2\H\^bcs) 

= {QiQ2\Ha + Hc\VLBcs)- (33) 

The three contributions to the gap equation are illus- 
trated in Fig. 5. In the next section we will write down 
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6m 





the explicit form of the gap equation and discuss the nu- 
merical solution. 



FIG. 5: The three contributions to Eq. (p2|). Smo is de- 
termined by the kinetic term, 5mc by 14// and SrUg is the 
contribution of the QQG cluster. 



III. QUARK MASS GAP 

From translational, rotational and global color invari- 
ance of the vacuum it follows that for each quark flavor, 



rfk 
(2^ 



.5('"?)(|k|)fet(k,A„z,, 



CT-k 



(34) 



The chiral angle is given by (c/. Eqs. ([191)), 

25('?'?) _ 2S'«9(|k|) 



tank = tan0g(|k|) = 



1-|S'99|2 1- (S''?9(|k|))2 



(35) 

To evaluate the matrix elements in Eq. (|21]) the Hamilto- 
nian needs to be expressed in terms of the quasiparticle 
operators. This can simply be done by noticing that in 
the quasiparticle basis the field operators become, 

^.(x)= E ^ 

A=±l/2 



(27r)3 



[C/(k,A)i3(k,A,z) 
+V{-k, X)D\-k,X,i)] 



(36) 



where the quasiparticle spinors U and V are given by 
1 r {E + M)x{X) 



C/(k,A) = 

v{-KX) 



y/2E{E + M) V ■ kx(A) 

1 f -ylCT-kx(A) 



^2E{E + M) \iE + M)x(A) 



(37) 



with E = £:(|k|), M = Af(|k|) = £;sink, A = A{\k\) = 
iJcoSk. Here we have introduced an arbitrary function 
£'(k) to make the expression for the single quasiparticle 
wave functions analogous to those of free particles, but 
it is clear that U and V do not depend on E but only on 
the chiral angle. 

Similarly for the gluon fields we have, 
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E 

A, a 



^2 '''^l'^2 



and in terms of the quasi-gluon operators the fields are 
given by, 

A'^(x) = E Z ,,!^,, J^ [a{k,X,a)e{k,X) 



\=±i 



-fa^(-k,A,a)e*(-k,A)] 



with 



and 

tanhk = tanh 



|k|(coshk + sinhk), 

25'(99) 2599(|k|) 



(39) 



(40) 



1 



1 599 I 



1 



(41) 

Truncating S at the QQG level leads to uncoupled gluon 
and quark gap equations. The gluon gap equation was 
studied in Ref. The gluon gap function ti^dkl) was 

determined by the matrix element of the Coulomb oper- 
ator in the BCS vacuum, which in turn was selfconsis- 
tently determined by the gluon mass gap. It was found 
that a good analytical approximation to, Ve//(x— y) (c./. 
Eq. (E4) ) is, in momentum space, given by. 



eff[ 



/(k)rf^(k) 

k2 



(42) 



where d(k) is the expectation value of the Faddeev-Popov 
operator and it is approximately given by. 



d(k) 



3.5 



f log(2.41 
Vlog(l-4 



0.48 



41 + |k|Vm?: 



and 



5(99)(|k|)at(k,A,a)at(-k,A,a). (38) /(|k|) ^ 



(w) 



0.97 



f log(l 
\^log(0.82+| 



82) 



for |k| < TOg 
for Ikl > TO„ 



for Ikl < 



for |k| > 



(43) 



(44) 



originates from renormalizing the composite Coulomb 
kernel. The gluon mass, rUg arises from dimensional 
transmutation and can be fixed by the string tension. 
The result of the fit to lattice data, yields nig ^ 1.6/ro ~ 
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600 MeV and is show in Fig. 6. The gluon gap function 
t^dkl) is weh approximated by, 



and 



u;(|k|) 



TUg for |k| < nig 
Ikl for k| > TO„ 



(45) 



The first two terms in Eq. ( |32| ) are then given by 

(5too = 5mo{\q\) = |q| sinq, (46) 



Smc = 5mc{\(i\) 



Cf 
2 



dk 

(2^ 



Ve/zdk - q|) sinkCOSq -k • qsinqCOSk 



(47) 



The new contribution to the gap arising from the QQG 
cluster contains matrix elements of H evaluated between 
the BCS vacuum and a tree particle, QQG state or be- 
tween QQG and QQ states. Only Vqg and He con- 
tribute to those and they are of order 0{g ^ ((i(|k|))) and 
0{g^ ^ (d^(|k|))) respectively. As discussed in Ref. Q 
the later is of a type of a vertex correction and is expected 
to be a small 0(10—20%) correction to an 0{g) contribu- 



J 



tion from Vqg. Therefore we will not further included it 
here (this is also consistent with the ring-rainbow approx- 
imation to Veff). The final expression for SrUg also re- 
quires QQG wave functions i.e the eigenstates of Hq+Hc 
projected onto the QQG states. In this work we do not 
attempt to solve this eigenvalue problem instead we will 
approximate the sum over 3-particle intermediate states 

by, 



El 



En 



E 



\QQG) 



QQG,{EQQG-Enscs>^i^ 
I 



E, 



QQG 



En, 



-{QQGl (48) 



i.e. we approximate the sum over the complete set of 
eigenstates by a single state with energy smaller then 
some factorization scale, Ap, E^, — Eq^^^^ < A^, and a 
perturbative continuum of states with energy grater then 
Kp. The scale Kp should roughly equal the energy where, 
due to string breaking the linear confining potential satu- 
rates. For the first excited hybrid potential Af > 1.5GeV 
which corresponds to the distance between color sources, 
r = 1.15 fm ||l^. Thus we expect that the size of the 
momentum space wave function, {QQG\^) should be of 
the order /? ~ 1/r = 0.2GeV. As for the spin-orbital 
momentum dependence of the QQG wave function we 
shall assume that it corresponds to low values of the or- 
bital angular momenta which are consistent with those 
of the low lying gluonic excitations in the presence of QQ 
sources. Lattice computation of the QQ adiabatic poten- 
tials arising from excited gluon configurations indicate 
that the so called n„ potential has lower energy then the 
Ilg potential [|l2| . These two correspond to gluon config- 
uration with J^'~^ — 1"' and 1 respectively which is 
also consistent with the bag model representation of glu- 
onic excitations [pSl. The QQG wave function coupled 



with the J^'^ — 1^ gluon quantum numbers would also 
have the QQ pair with the same J^'~^ = 1^ quantum 
numbers (to give the overall J^*^ — 0++ of the vacuum) 
and would be given by 



[QQG]' 



X [(Lg - 1) X (5g = 1)]' 



(49) 



It is easy to check, however, that since Vqg is spin depen- 
dent this wave function has vanishing overlap with the 
Vqg\^) state. The other possibility is to take J^'-' =1^^ 
configurations for both the glue and the quark-antiquark 
which give, 

[QQG]' = [[(Lqq = 0)x(5qq = 1)]^ 

X [(Lg = 0)x (5g = 1)]']". (50) 

and take the spin-orbit wave function in the form of. 
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FIG. 6: Comparison of the lattice results for the ground state potential between two static 3 and 3 sources and the fit to Ve// 
given by Eq. (^) (solid line), and Eq. (^) (dashed line). Lattice data (boxes) come from Ref. ||l^ 



(Q(kQ, Aq, ig), 0(kQ, Aq, iq), G(kG, Ag, = (2^)35(kQ+kQ+kG)C/^(kQ, AQ)a-e(kG, AG)l^(kQ, AQ)vI/(kQ, k^, kG). 

I m 

The color part of the wave function is given by T°'^ , The expression for Srug is then given by, 
and for the orbital wave function we will take a gaussian 

ansatz, Snig = Smg{\q\) = Srug^soft + Snig^hard, (53) 

M'(kQ,kQ,kG)=exp(-(k^+k|+k2.)//32). (52) with 



6m 



g.soft 



Cf 



dk d(|k-q|)vl/(k,q,k-q)/|*| 



- En J (27r)3 ^cj(|q-k|) 
dp dl d(|l-p|)**(p,l,l-p)/|*| 



[SkCq - /(k,q)cuSq] 



(27r)3 (27r)3 



V^(ll-Pl) 



[l + SpSl+/(p, l)CpCl], 



(54) 



Sm 



g.hard 



-Cf 



dk d2(|k„q|) (I_vl/A^(k^q,k_q)) 



(27r)3 |q-k| 



|k-q| 



SkCq - /(k, q)CkSq] . 



(55) 



Here 

(|kp + |qp)k.q-|k||q|(l + (k.q)^) 
/(k, q) = — — — , (56) 



UV behavior and renormalization 



dp dl 

(27r)3 (27r)3 



(k-q)2 



|*(p, 1, l-p)|2 [1 + SpSl + /(p, l)CpCl 



and vE'Af. given by Eq. (|5|) with /3 ^ A^;- so that 1 - *Af 
cuts off hard QQG contribution for energies below Ap. 



Before analyzing the full gap equation and in particu- 
lar the effects of Srug, we shall first discuss the IR and 
UV behavior in the BCS approximations. The BCS ap- 
proximation to the chiral gap has been studied earlier 
for various model approximations to Vef / . Most of them 
(57) use an effective potential which is regular at at the ori- 
gin, e.g a pure linear potential Ve//(r) = 6r |^ or a 
harmonic oscillator, Veff{r) = kr^ B. For such poten- 
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tials the gap equation is finite in the high momentum 
Umit and no renormaUzation is required. This is not 
the case if potential has the Coulomb component with 
Veff{r — *■ 0) ~ a/r and a being either a constant or 
a running coupling a — > a(r) ~ l/ln(l/r). The BCS 
quark gap for potentials with the Coulomb tail was stud- 
ied in Ref. ||, Q and Ref. joj. The gap equation used 
in Ref. 1^ would be identical to one used here, if 5mg 
was set to zero {e.g. the BCS approximation). Instead, 
in Ref. ||^ an energy-independent interaction motivated 
by a transverse gluon exchange was added. In Ref. ||] 
it was argued that, in the chiral limit, the resulting gap 
equation, could be renormalized by introducing a single 
counterterm representing the wave function renormaliza- 
tion. Starting from the Coulomb gauge Hamiltonian this 



would arise if the free quark kinetic energy term was re- 
placed by a renormalized one. 



dxi/)(x) [—ia ■ V?/;(x)] 

Z(A) / dx?/;(x) [-ia ■ Vi/'(x)]a (58) 



The explicit, UV cutoff- A dependence regularizing the 
kinetic operator can be introduced, for example by field 
smearing, however, the regularization procedure becomes 
irrelevant once the resulting gap equation is renormal- 
ized. The unrenormalized BCS, gap equation (without 
effects from transverse gluons) is then given by. 



J 



Z(A)m(|q|) = 



Cf 



— l-e//(k-q) ^ 



Cf 



i2n) 



^k2 + m2(k) 2 J (27r) 
I 



dk , Ikl m(q) 
— K//(k - q)k . qU W_ 



q| ^k2 +m2(k)^ 



(59) 



where we have defined the constituent mass, rn(|k|) by solution, m(|k|) at a specific value of |q| — |qo|. This 

r7i(|k|) = |k| sink. The renormalized equation is obtained leads to a (A and qo-independent), renormalized gap 

by a single subtraction i.e. by fixing the A-independent equation, 

I 



m(|q|) lim [/™(|qo|. A) - /z(|qo|, A)] = TO(qo) lim [/„,(|q|. A) - /z(|q|, A)] , 



(60) 



/m(|q|,A) 



Cf_ 
2 



dk 

(2^ 



-K//(k-q) 



TO(k) 



^k2 + m2(k) ■ 



^z(|q|,A)^ 



Cf 



V 



dk , |k 
— K//(k-q)k.q^ 



(2^) 



m(|q|) 



q| ^ 



i611 



Whenever possible we will also use the notation /(|q|) = 
/(|q|, oo). We will now show that, this equation does not 
have a well behaved, nontrivial solution vanishing asymp- 
totically in the large momentum limit, as it was assumed, 
for example in Ref. Before we do that first we need 
to take care of the possible IR divergences which appear 
in the integrals when k ^ q. In this limit Ve/ / (k — q) is 
highly divergent, reflecting the long range nature of the 
confining interaction, e.g. Ve//(k — q) cx l/(k— q)* for 
the linear potential. The gap equation, however, is finite 
due to cancellation of the numerators between and 
Iz- To make individual integrals well behaved we can 
split the lU and UV parts of Vejf defining, 

y^«(k,M)^0(M-|k|)Fe//(|k|), 

T/^^(k,M)^0(|k|-M)l^,//(|k|), (62) 

since y^^(k, M) + F^^(k, M) = V;//(k) and gap equa- 
tion is independent on the parameter M and we will not 



write it explicitly. The gap equation becomes, 
1 



m(|q|) 



^ + s(|q|) 



|q|) 



where 

A = A(|qo|) 



^r(|q|)] , 

(63) 



[/4^(|qo|)-/F(|qo|)+4r(|qo| 



TOllqol 



and 



i3(|q|) = S(|q|,|qo|) = 



jUV 



(|qol) 



"^(|q|) n^(lqol) 



(64) 



(65) 



For given qo , >1 is a constant and S is a function of q, and 
both, A and B are well defined. In A the IR divergences 
cancel between /4^(|qo|) and /^^(|qo|), and J,^^ is finite 
if m(|q|) — > as q ^ oo. In B each term is IR finite 
and the UV divergences cancel between the two terms in 
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Eq. (|6|). It is easy to show that for |q| >> M, |qo| the 
function i3(|q|) , behaves as, 



-Cf— log q 
37r 



(66) 



for Ve//(|k|) 47ra/|kp as |k| oo. If a is replaced 
by a running coupling then |_B(|q|)| grows with |q| like 
log log q^. From Eq. (|6^) it thus follows that for some 
|q|/|qo| >> I, A + -B(|q|) changes sign and therefore 
the equation is undefined. This also remains true if an 
additional transverse potential is added as done in 
In this case the the argument of the integrals defining 
hmction the -B(|q|) becomes, 



k.qgy^^(k-q) 



M 

|q| 



l^^^(k-q)k-q 

+2/(k,q)F^^^(k-q)] 

(67) 



For Vrik) = 47ra/(k2 + const) (as used in Ref. ||) the 
additional transverse potential does not contribute to the 
logq^ (or log log q^) behavior of i?(|q|). In our case there 
would be a similar contribution arising from the hard part 
of the gluon exchange given by Sirig^hard- At large |q| it 
adds, a positive , +CFYfjlogq^ contribution to i3(|q|) 
and therefore does not cause problems on its own but 
at the same time does not eliminate the singularity from 
the Coulomb potential since the net effect is such that 
_B(|q|) —oo as |q| oo. 

The problems with the renormalized gap equation for 
the Coulomb potential is illustrated in Figs. 7 and 8. In 
this test case we simply take 



Veffik) = 



V^^^(k) + T/^^(k) 

1 Snb Ana 
C^^^ k2 log(k2/m2 



2)" 



(68) 



For the string tension, b = 0.24 GeV^, a = 0.1, and n — 
this gives a good fit to the lattice data as shown in Fig. 6. 




q [GeV] 

FIG. 7: Solution of the gap equation for a potential given by 
Eq. (H), with a = 1, and n = 0, 1, 1/2, 3/2. The three upper 
lines correspond to n = 0, the next three to n = 1, n = 1/2 
and n — 3/2 respectively 
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In Fig. 8 we show the function B{q) calculated from a 
numerical solution to the gap equation, Eq. (|63|). Since 
Veff is already given by a sum of two terms, one dom- 
inating in the IR and the other in the UV the compo- 
nents V^^ and V^^ can be defined as the linear and 
the Coulomb piece respectively. If n > 1 there is no 
renormalization required and the gap equation is given by 
Eq. ( |59| ) with Z — 1. The n — 3/2 case corresponds to an 
approximate analytical solution for Vef / discussed in |l^ 
and is also close to the exact, numerical solution given 
by Eqs. ( |42| ). The solution of the gap equation, m{q) for 
n = 3/2 is shown in Fig. 7 by the lowest line (boxes). The 
function -B(|q|) corresponding to this case is shown in 
Fig. 8 by the upper solid line (at large |q|), which asymp- 
totically approaches _B(|q|) — > const— 2CFQ!/(37rlog(q^)) 
as |q| oo. We then take this solution to set the value 



FIG. 8: Function -B(|q|) calculated for the potential as in 
Fig. 7. The lines (from top to bottom at high |q|) correspond 
to n = 3/2, n — 1/2, n = 1 and n — respectively 



of 7Ti(|qo| = nig — 600 MeV) and solve the renormalized 
gap equation, Eq. ( |63| ) for n = 1, n = 1/2 and n = 0. 
The asymptotic behavior at large |q| of i?(|q|) for these 
three cases is given by, 

( -CFi^\oglog{\q\^), n = l 
S(|q|) ^ <^ -2Cf|; W/'dql), n = 1/2 (69) 
[ -CF^log(|q|), n = 

The corresponding solutions to the gap equation are 
shown by the five upper lines in Fig. 7. The highest 
three correspond to n = case and their splitting in- 
dicates that the numerical procedure has not converged 
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into a unique solution. These three solutions correspond 
to three different cut-offs on the maximum momentum, 
|q|maa; = lOrUg, lOOrUg and lOOOmg. The other two lines 
correspond to solutions for n — 1/2 and n = 1 respec- 
tively. In these two cases the same three values for the 
momentum cutoffs were used and apparently in both 
cases a cutoff independent solution has emerged. This is 
because for n = 1/2 and n = 1 \B{q)\ grows very slowly 
and in practice the zero of the denominator in Eq. (^) 
is not crossed. This test calculation was performed with 
unphysically large a — 1. For a < 0.5 numeral com- 
putations, which always have a build in a finite upper 
momentum cutoff converge for |q|maa; as large as lO^irig. 

It is clear that the problematic UV contributions orig- 
inate from need for wave function renormalization. This 
problem has been resolved in Ref. H using an effective 
Hamiltonian with perturbative 0{g ) contributions cal- 
culated via a similarity transformation |p^ . In that ap- 
proach, in addition to the Coulomb and transverse gluon 
contributions to the gap equation, Smc and Srug^hard, 
there was also a modification of the single particle ki- 
netic. The additional contribution to the gap equation 
via 6mo cancels the log|q| term from -B(|q|) and results 
in a well defined equation. The disadvantage of that ap- 
proach however, is that it is restricted to the free rather 
then BCS basis and so far it has not being generalized 
beyond perturbation theory. 

The resummation of the leading UV contribution 
to the Faddeev-Popov operator and the Coulomb ker- 
nel has the effect of softening the UV behavior (c./. 
Eqs. (HI), (HI)) and at the BCS level leads to a finite gap 
equation without need for any additional, e.g. wave func- 
tion renormalization counter-terms. Furthermore the 
method enables to include effects of transverse glu- 
ons with a well defined energy dependence. The large 
momentum contribution from the QQG cluster, Smg^hard 
still requires renormalization through the wave function 
counterterm. As discussed above since it leads to -B(|q|) 
which is positive at large |q| the renormalized gap equa- 
tion is well behaved. 



0.15 



B. 



Numerical results 



We will now discuss the numerical results. These are 
summarized Fig. 9. The BCS, potential contribution to 
the gap equation, Smc is split into an IR and UV parts 
by setting /(k) = = d{k) = for |k| > mg and |k| < 
rUg respectively. For the pure IR potential (dashed line 
in Fig. 9) the gap function, m{q) is below 100 MeV for 
low |q| < TUg and it vanishes rapidly (as l/|q|"') at high 
q. The addition of the UV component of the potential, 
i.e. the Coulomb tail with the l/\og{q)^,n > 1 UV 
suppression, does not change much the low momentum 
behavior of m(|q|). It actually decreases m(0) to about 
75 MeV, by it increases the high momentum tail, overall 
leading to no change in the {QQ) condensate which stays 
at about -(111 MeV)3. 




0.05 - 



10" 10' 
q [GeV] 

FIG. 9: Solution to gap equation. The dashed line is a 
solution with the IR part of the potential (linear potential) 
only. The lowest solid line comes from a solution using the full 
static potential, Smc. The next higher solid line includes the 
static potential and the hard gluon, Smgj^ard, contribution 
from the QQG sector calculated for Af ~ 1.5 GeV. The 
shaded region corresponds to the full solution with 0.1 GeV < 
P < 0.4 GeV, and 1 GeV < Af < 2 GeV 



The Smg contribution depends on, f3 which sets the 
size of the soft wave function, Ap which divides be- 
tween the soft and hard one-gluon intermediate states, 
and Eq, — which determines the energy of the low 
gluonic excitations. As discussed earlier it is reasonable 
to set P ^ 0.2_GeV and Af ~ 1.5 GeV. As for the energy 
of the soft QQG state we take, Eq, — Eq_ = mg which we 
expect to be close to the lower bound and would therefore 
give the upper limit on the QQG contribution. The effect 
of the hard one-gluon-exchange contribution, defined by 
Smg^hard, is to iucrcasc m{q) yielding m{0) ^ 80 MeV 
and enhancing the condensate, {QQ) = —(150 MeV)'^. 
The solution to the gap equation including Smc and 
Smg^hard is shown by the second to lowest solid line in 
Fig. 9. 

As mentioned above, the gap equation with Smg, hard 
requires renormalization and we have simply set Z — I at 
A = |qmaa;|. Usiug any of the three values of |qmaa:| given 
previously no effect on the solution could be observed. 
This is analogous to the test case discussed earlier. 

The full effect of the QQG sector, including the soft 
contribution, parametrized by Smg^soft with the factor- 
ization scale /3 in the range between 0.1 GeV and 0.4 GeV 
and Kp = 1 — 2 GeV is shown by the shaded region. The 
lower limit corresponds to /3 = 0.1 GeV and A^ = 2 GeV 
and yields {QQ) = -(140 MeV)^^ for the upper hmit 
(3 = 0.4 GeV, Af = 1 GeV and (QQ) = -(190 MeV)^. 
The addition of the soft gluon intermediate state brings 
both, the constituent quark mass and the condensate sig- 
nificantly closer to phenomenologically acceptable values. 

An alternative simple parameterization of the soft 
gluon contribution would be to replace it by an effec- 
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tive local operator, by expanding Snig^soft in powers of 
|k|/TOg. The lowest dimension operator has the structure 

= [^A^(x)aV(x)(5T,.,(Vx)^^(x)a-'"^(x)]^ 
a 

.(70) 

with C being a dimensional constant, l/w^ scale arising 
from the product of — Eq^^^ and w(0) = nig, and the 
operator being related by the factorization scale Ap. The 
appearance of these to scales is quite natural. Since the 
operator arises through elimination of part of the Fock 
space the overall scale is given by the excitation energy of 
the eliminated sectors and the momentum cutoff comes 
from the spatial extend of the excited state wave func- 
tion. Such a simple, local approximation of the soft QQG 
exchange was considered previously in Ref. [20j where it 
was shown that such an operator was indeed relevant to 
chiral symmetry braking effects, e.g. the condensate and 
the TT — p mass splitting. 



Coulomb interaction. However, the extent of chiral sym- 
metry breaking generated this way, as measured by the 
scalar quark density or as compared to the phenomeno- 
logical constituent quark model is to small. We have 
shown that the naive inclusion of the short range part of 
the QQ potential leads to instabilities in the quark gap 
equations which cannot be renormalized away. In con- 
trast a systematical re-summation of the leading IR and 
UV corrections to to the bare Coulomb kernel leads to an 
effective interactions which is consistent with the varia- 
tional treatment and the gap equation. Using the linked 
cluster expansion we have estimated the role of three- 
particle, QQG correlations in the vacuum, and shown 
that they are indeed important, in particular their low 
momentum components. Even though we have not used 
the exact solution describing the soft QQG state our re- 
sult are expected to be close to the upper bound for the 
non-BCS contribution to the chiral condensate and are 
consistent with previous studies. 



IV. SUMMARY 

Dynamical breaking of chiral symmetry is the funda- 
mental property of QCD which leads to the constituent 
representation. In the canonical, Hamiltonian based 
formulation it arises via the BCS-like pairing between 
light quark and antiquarks mediated by the attractive 
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